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Abstract 

, ^ The variational approach, used by Feynman in the study of the polaron problem, is generalized 

to treat a system of N non-relativistic particles interacting with scalar and vector mesons. After 
integrating out the meson fields in the path integral formulation we perform a variational calculation 
for the elfective N-body two-time action to obtain the energy and the mass of the system. The 
interplay of self energies and exchange terms in attaining binding and saturation is examined. We 
estimate the size of the particles in the medium and give the mean number of mesons as a function 

Cn of N. 
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1 Introduction 



The description of a strongly interacting many body system starting from the underlying microscopic 
theory is a fundamental problem in theoretical physics. It is well known that we need a relativistic 
field theory to describe the fundamental interactions of such an N-body system. However solving for 
the N-body bound state properties in an exact field theoretic approach is still an unsolved problem. In 
traditional nuclear physics the approach has been to use in the nonrelativistic many body Schrodinger 
equation static or near static nucleon-nucleon potentials, deduced from nucleon-nucleon scattering. 
The solution of the resulting many body Schrodinger equation in most cases involves a number of 
approximations such as the summation of a selected class of diagrams as in the Brueckner approach. 
Application of Monte-Carlo techniques, which proved to be extremely useful in the study of single 
particle properties such as masses of mesons and baryons [1], remains a major problem in the case 
of many fermions, due to antisymmetry. It is only recently that variational Green's function Monte- 
Carlo techniques have began to yield a more complete solution to the non-relativistic Schrodinger 
equation in the case of selected nuclei [2]. Despite the success of the traditional approach we know 
that it does not present the full picture. We know that the interactions arise from meson exchange 
leading to non-static potentials, vertex corrections and self energies and that mesons can have a real 
existence in the nucleus. Therefore one would like to study the many-body system with dynamical 
mesons present. This is the main motivation of the present investigation. 

In a previous work [3] the dressing of a single nucleon by mesons was examined using as underlying 
theory the Walecka model [4]. In this model the main features of the nucleon-nucleon interaction arise 
from the exchange of scalar and vector mesons. The method used was based on Feynman's variational 
approach [5] originally applied to the Frohlich Hamiltonian [6] describing an optical polaron. This 
approach has been shown by a numerical calculation [7] to be the best semi- analytical method valid for 
both weak and strong electron-phonon coupling. The success of the method is based on the fact that 
the retardation effects of the interaction in the original system are taken into account by simulating the 
system with a retarded harmonic oscillator with properly chosen coupling and retardation constants. 
In order to apply this approach to the Walecka model we treat the nucleons, being heavier than 
the mesons, nonrelativistically and keep only terms of order 1/M where M is the nucleon mass. To 
this order there is also a spin-dependent part which we neglect for technical reasons. Within the 
nonrelativistic approximation the Walecka model resembles closely the Frohlich Hamiltonian. In the 
path integral formulation the mesonic degrees of freedom can be integrated out in an exact manner 
giving rise to a non local (two-time) effective N-body action. 

The properties of the N-body system are evaluated using Feynman's variational ansatz generalized 
for the N-body action. In this way we sum up self energy and exchanged diagrams to all orders, albeit 
in an approximate way, with only the first order perturbative results reproduced exactly for the energy 
and effective mass. In addition to the binding energy and mass of the system we calculate the size of 
the dressed constituent in the medium and the excess number of scalar and vector mesons in the bound 
N-body state, both of which are interesting quantities in connection with the various explanations of 
the EMC effect [8]. The results obtained throughout this work should only be taken as a qualitative 
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description of the corresponding quantities of real nuclei. The reason for this is that, in addition 
to the nonrelativistic approximation and to neglecting the spin dependence, we have not imposed 
antisymmetry. As we will explain in section 3 imposing antisymmetry destroys the quadratic form of 
the trial action and one can no longer perform the functional integrations analytically, a prerequisite 
for the present approach. Thus the strong repulsive short range interaction is the only factor that 
prevents two particles overlapping with each other. We expect the effective "hard" core to be softer in 
comparison to the real nuclear case. Nevertheless the bound states approximate a saturating system 
for special values of the coupling constants and we will still refer to the constituents as nucleons. 

The paper is organized as follows: In section 2 we derive the many body action for the Walecka 
model after elimination of the mesonic degrees of freedom and discuss the static two-body potential. 
In section 3 we calculate the energy and the effective mass of the multiparticle state. We derive also 
expressions for the mean number of scalar and vector mesons and for the radius of the bound state. 
In section 4 we present our numerical results and discuss the stability of the N-particle system. The 
details of the evaluation of the various observables can be found in the appendices. 



2 Effective many-body action 

The interaction part of the Langrangian in the Walecka model consists of two contributions: i) an 
attractive interaction between the nucleons (tp) due to the exchange of scalar meson fields (cj)) which 
is responsible for the binding of the nucleons and, ii) a repulsive short range term which is described 
through the exchange of vector meson fields (F^) and is responsible for the stability of the multi- 
nucleon bound state. In terms of the nucleon and meson fields the Langrangian can be written as 

£ = £o(V', V^) + gsi^i^4> - gvh^^W (2.1) 
where the free part is given by 

£o(V', V^) = i,{i-i^d^ - M)i, + ^(9^<^9^<^ - m|<^2^ - \f^.F^'' + \mlV^V^ . (2.2) 

To obtain the effective action for the interacting nucleons we have to perform the path integral 
over the meson fields. There are no interaction terms between the scalar and the vector mesons in 
the above Langrangian and therefore the integration can be easily done leading to an effective action 
with two interaction terms which correspond to the scalar and vector contribution respectively. The 
exact form obtained is given by 

S = So + ^ J d'^xd'^x'[-glps{x)Gs{x - x')ps{x') + gvJf,{x)Gy{x - x')Jj,{x')] (2.3) 

where 5*0 is the kinetic term, ps{x) = tp(x)tjj(x) is the scalar density and = ipixYj ^ipix) is the 
baryon current. Gy and Gs are the propagators of the vector and scalar mesons respectively. The 
scalar propagator Gs is given by 

„ /• d'^k ie-^k.(x-x') 

Gs(x-x')= - — -— 5 (2.4) 

^ ^ J (27r)4 A;2 - ml + le ^ ^ 
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and the corresponding expression for the massive vector field propagator, G'^ , reads 

It is well known [9] that, due to the k^k'^ term, the vector meson propagator develops a singular ultra- 
violet behaviour. However the baryonic current is conserved and we can therefore neglect this term, 
since, in every physical quantity, it always occurs contracted with a current of vanishing divergence 
[4]. 

In order to stay as close as possible to the Feynman ansatz the main assumption that we make is 
to consider the nucleons much heavier than the mesons so that we can treat them non-relativistically. 
We realize the non-relativistic approximation by going over to the particle picture for the nucleons 
keeping terms of (9(1/M). The scalar density is then given by 

N 

ps{x) ~ ^l}\x)^l}{x) = ^ ^(x - x^) (2.6) 

8 = 1 

and the vector current [10] by 

dx ^ 

.J^'{x) = ^.J'^{x) , /f (a;) = -j^^3(x- x,(i)) + spin current , a;? = a;^ = ... = a;^ = i (2.7) 



1=1 



dt 



The above current is conserved, i.e. it fulfills the continuity equation. The non-relativistic approx- 
imation is incorporated in the use of a common time variable for all nucleon trajectories. The spin 
dependent part in eq. (2.7) will be neglected in what follows in order to avoid the complications of 
spin in a path integral [11]. Inserting the non-relativistic current into eq. (2.3) and performing the fco 
integration we get for the effective action 



N 



s. 



eff 



E 



d^k 



47r2 



Rs{k,t- t') Rv{k,t- t') 
ujvik) ^ i2y(A;,0)e*-(^'W-^^W) | 



f2.8) 



with 



R.ik, T)=f- exp i-io,ik) \t\ ) io.ik) = + m2 t = S,V. ( 



47r (jJiik) V 
We stress that up to this point the mesonic degrees of freedom are treated in an exact fully rela- 
tivistic manner. Considering the Walecka model as a low energy effective theory we have introduced 
form factors Fi(k) at the meson-nucleon vertices to cut off the high-momentum components. Here we 
use a standard monopole-type form factor [12] 

A? 



A2 + k^ ' 



i = s,v 



2.10 
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In deriving eq.(2.8) we have not respected the Pauli principle for the nucleons for purely technical 
reasons. To project onto the antisymmetric space within the present formulation one uses antisym- 
metric coordinate states [13] to obtain a determinant in the particle labels. The simplest way to see 
this is to consider a free gas of fermions of mass m. Then the infinitesimal evolution operator at the 
A;*^ time slice between antisymmetric coordinate states is given by 

/ \ N/2 

< xi(A;),...,XAr(A;)|e-^^|xi(A;- l),...,XAr(A;- 1) >= — DetM(k) e-f.}ZSMk)-Mk-i)? 

(2.11) 

where 

M^^^j.) ^ g-f [(x,(fc)-x,(fc-l)f-(x,(fc)-x,(fc-l)f ] 

and e is the infinitesimal time step. Therefore including antisymmetry on the level of the trial action 
introduces determinantal dependences on the particle trajectories and we are no longer able to perform 
the functional integrals analytically destroying one of the criteria required by our method. There are 
attempts to model the Pauli principle classically [14] and one could envisage applying those ideas 
within the present approach. Progress in this direction will be reported elsewhere. 
This effective action has several important features: 

i) Due to the interaction with relativistic mesons the k-dependence of the action is more compli- 
cated as compared to the polaron action where integration over the Ijk form factor gives rise 
to the retardated Coulomb potential (more on this point in appendix A). 

ii) The scalar mesons give rise to a two time effective action which closely resembles the polaron 
action for N = 1. 

iii) The vector mesons generate in addition a transverse velocity dependent contribution to the two 
time effective action which is absent in the Frohlich Hamiltonian. For N > 1 they also produce 
an equal-time exchanged term. 

The static two-body potential can be easily obtained from ^efr in the following manner: We make 
the transformation 

T = —— T = t-t' (2.12) 
2 ^ ^ 

and we consider the limit r ^ so that Xj(i) — Xj(i') = Xj(r + f) ~ ~ 2 ) ~ ^«(^) ~ ^i(^)- 

the absence of the form factors and using the integration formulas of the modified Bessel functions 

given in appendix A we can perform the r integration. We find the familiar sum of a repulsive and 

an attractive Yukawa potential, typical for vector and scalar massive meson exchange 

V(r) = ^ — ^ , , (2.13) 

^ ^ Att |r| Att |r| ^ ^ 

where r = Xj- — Xj. To appreciate the modifications introduced by the presence of the form factors we 
give the static potential for the monopole-type form factors used. It reads 



{(A-l-miy |r| 2(A2,-m2,) 
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where 

a, = ^ i = S,V. (2.15) 
47r 

If the cut off A5 and Ay are large compared to the meson masses ms and my then we obtain the 
result of (2.13). From fits to nuclear properties the values obtained for the meson masses are 

m5 = 520 MeV my = 780 MeV . (2.16) 

Since we have treated the Walecka model as a low energy effective theory for the nucleons the cut offs 
should be of the order of the nucleon mass in order to eliminate large momenta. This requirement 
in combination with the fact that the cut offs should be larger than the meson masses makes the 
non-relativistic approximation not very well fulfilled in this model. It would be more appropriate to 
study the pion-nucleon system within this approach if it were not for the complications of spin and 
isospin. We show the form of this potential for various values of the cutoff parameters A5, Ay in fig.l. 
As it can be seen, for this allowed range, the effects of the cut off are rather strong altering the height 
of the repulsive core and the depth of the minimum. 



3 Trial action and observables 

A good trial action must incorporate the essential physics contained in the effective theory as given by 
(2.8). At the same time it must be simple enough so that we can perform analytically the functional 
integrals necessary to calculate our observables. This gives us the following conditions for the trial 
action: 

i) It must be translationally invariant. 

ii) It must include the effects of retardation. 

iii) It can only be at the most a quadratic function of the particle trajectories. 

The last criterion will not be fulfilled if antisymmetry is included. To construct a suitable trial action 
we follow Feynman's approach and represent the mesons by fictitious particles harmonically coupled 
to the nucleons. To take into account the self energies we attach a fictitious particle by a spring to 
each nucleon. We will refer to this particle as the y— particle. The exchange of mesons is modeled by 
attaching all nucleons to an additional fictitious particle, the Z— particle. This system is schematically 
shown in fig. 2 and corresponds to the Langragian 



^ '"M.o m„ . 9 ii' „ ,9 K, 

yx^ + + - _(z - X,)' - -(y. - X, 



E 



=1 



(3.17) 
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Integrating out the fictitious particles in the Feynman path integral we obtain after rotation to Eu- 
clidean time and taking the (3 ^ oo limit 

+ dt'{^,{t)-^,{f)fe-^\'-'\\ , (3.18) 

with 



Since the above action is quadratic in Xj- all the path integrals can be done analytically. However 
the resulting expressions for the various observables are rather involved and are given in appendix 
B. Instead we consider here a special case of the above trial action derived from the Lagrangian of 

(3.17) by setting k = 0; in other words we model both exchange and self-energy contributions by a 
single fictitious particle. It corresponds to setting C'd = in (3.18) so that exchange and direct terms 
have the same strength and retardation reducing the number of variational parameters to two. It is 
the simplest extension of Feynman's trial action for the many body system and we will refer to it in 
what follows as trial action 1. This trial action is used mainly as an illustration of the main steps 
of our calculation. In order to demonstrate that it is important to have different strengths for the 
diagonal and non-diagonal terms in the trial action irrespective of the retardation we consider the 
special case of allowing different C and C'd but still keeping the same retardation. It corresponds to 
a model Lagrangian of the system shown in fig. 2 but with K/niy = NK/Mz and it emerges from 

(3.18) if we set w = W£,- We will refer to this 3-parameter action as trial action 2. It allows us to 
investigate the relative strength of diagonal and exchange terms of the trial action. This separation is 
lost in the four parameter trial action (trial action 3) given by (3.18) since, in addition to the purely 
diagonal terms with retardation W£,, there are diagonal terms having retardation factor and strength 
the same as that of the exchange. A way to recover this separation is to exclude diagonal terms in 
the sum of the exchange terms in Sp, leading to trial action 4. This will be useful when we discuss 
the dressing of the constituent which is the main reason of introducing trial action 4. In contrast to 
the rest of the trial actions it can no longer be derived from an underlying model Lagrangian. Trial 
action 4 completes the variational ansatze that we considered in this work. 

A variational bound on the energy is obtained on the level of an action by using Jensen's inequality, 

^ ; , — 5 =< e-^ > > e-<^> . (3.20) 

This gives for the energy of the multi-particle system the upper limit 

En<Ef+ lim ^ ^-^-J'f > (3.21) 

where Ep is the energy corresponding to the trial action Sp and the expectation value is taken with 
respect to the trial action Sp- 
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In order to obtain the effective mass Mes(N) of tlie iV-particIe system we calculate tlie contribution 
to tlie energy by giving tlie center of mass mean velocity Ucm- Tlie trajectories of the nucleons now 
satisfy the boundary condition Xj(/3) = Xj(0) + Ucm/3- Using the calculated energy _Ejv(ucm) we obtain 
the effective mass of the system from 



Meff(iV) = En{u. 



cmi|Mcm = 



f3.22l 



The path integrals involved in the computation of the expectation value in eq. (3.21) can be done 
easier using a Fourier representation for the particle paths with the path of the ith particle given by 



Mt) = + E cos(27rn-) + ^ sin(27rn-) . 

^ n=l P n=l P 



f3.23l 



As long as the particle trajectories occur at most quadratically in the exponent of the path integral 
the stationary phase method is exact and easy to apply. The details of the calculation can be found in 
appendix B. There we also give the results for the energy, the effective mass and the other observables 
using trial action 3 ( and consequently 2) and 4. Because of the complexity of these expressions we 
will discuss here the results obtained using trial action 1. For the energy we find 

Ej^ = Ef - - j dr j dkk^ y^Rs{k, r) i^N - l)e-^^^^^^ + e'^^^o^^^j + NRv(k, r) 



(N-l) 



m 
T2 



V 



c^e '"'^ - c V e 



MNv 



,2 

g 2M 



c^e-"^ + {N - l)c v e 



1 



dk N(N - 1) Rvik,0) ujy e"^ + Co 



where we have used the notation 



2M (''") 



e 2M f^D 



(3.24) 



UCN 
Mw ' 



2 2,2 
V = C + VJ 



(3.25) 



Ep is the energy corresponding to the trial action 1 and is given by 

3 - w f 



Ef 



+ (iV- l)c 



f3.26) 



The quantities l-i^{T) and l-ij){T) are given by 



1 v'^ - w'^ 



N 



1 I - 



N 



'l-e 



l-e- 



— r + - iV- 1 + e- 

f ^ c 



[N - 1) 



1 - e" 



(3.27) 



and Co, which is independent of the variational parameters, by 

A[ roo / 1 f. 



(3.28) 
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For the effective mass we find 



Mefr(iV) = Mn(i 



37r Jo 

2M 
( ^ 



I 2M 



Rs{k,T) 



i2y(A;,r) 



c f e 



MNv 



e 2M 



(iV - f)e"^ 



f3.29l 



A cfieck of our variationai caicuiation is tfie reduction of our expressions to tfiose obtained in 
perturbation tfieory wfien gv ^ and gs ^ wfiicfi in our formaiism means v w. Taking 
V = w -\- O(g^) in tfie above expressions tfie contribution from tlie non-diagonal part of tlie action 
vanislies and we obtain tlie sum of the single particle self-energies i.e. no binding which is the 
perturbative result. 

Other interesting parameters of the many particle system are the average interparticle distance 
and the mean number of mesons present in a bound state. The mean number of scalar and vector 
mesons is obtained by finding the expectation value of the corresponding number operator to zeroth 
order i.e. taking the expectation value with respect to the trial action. The method is the same as for 
the one particle case described in [3] and therefore we only give here the resulting expressions. We 
find for the mean number of scalar mesons 



< Ns >-- 



vr Jo 



dk I dr k TRs{k, r) 
Jo 



3.30 



and of vector mesons 



< Nv > 



vr Jo 



dk I dr k\Rv(k,T) \(N -1) e-^^ 



m 



V 



1 



Ll>v dfi^ir) 
2M dr 

2M dr 



1 d'n^T) 

M dr^ \ 
M dr"^ 



1 



2ml 



1 



f3.31^ 



The above formulae hold for a general quadratic trial action and for N = 1 reduce to those given 
in ref. [3]. The explicit expressions for fi'jj(T), I-i'^(t) and their derivatives depend of course on the 
specific trial action. 

The interparticle distance is defined by 

N 



r2 



1 



N(N 



— ^ < (x, - Xj 



> 



f3.32l 



«' j'=i 



and for a system showing saturation i.e. constant particle density the size of the system is related to 
the interparticle distance by 
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We evaluate the interparticle distance to zeroth order. The simplest way to do this is to use the wave 
function corresponding to the variational Hamiltonian. For the simplest trial Hamiltonian we have 
harmonic couplings of the nucleons to one fictitious particle. In this case the intrinsic wave function 
is given by 



-(/'(xi, ...,xjv, Z) = A/'exp 

where 



N n N 

-a^(Z-x.)^-|^ (x.-x,)^ 

i=l ^j' = l 



(3.34) 



a = ^VKM,,d /3 = l|^^:^-aj , (3.35) 

Mj-ed = m+m'z reduced mass and M is the normalisation constant. In addition to the interpar- 

ticle distance a quantity of interest is the size of the particle in the medium. If we were to probe a 
nucleon with a meson the nucleon - meson vertex will be dressed giving rise to a form factor. The size 
thus measured can be considered as the core radius of the nucleon. In the one particle case this defi- 
nition was shown [3] to be identical in zeroth order to the nucleon root mean square radius calculated 
using the intrinsic wave function of eq. (3.34). We use the same definition in this many particle case 
where we consider the fictitious and the ith particles as an independent subsystem. Relative to the 
center of mass of this subsystem we obtain for the root mean square radius or center of mass radius, 
I'cim of the nucleon in the medium 

The cloud radius, rdoud, is a measure of the meson cloud around the bare nucleon and in the one- 
particle case it is defined as the distance of a particle from the fictitious particle. The generalization 



1 A .9 3 ( V N -1 



r 



cloud 



(iV) = - g < (Z - x.)^ >=_!_ + ^=i= (3.37) 



for the N-particle system includes, in addition to the self energies, exchange contributions and it is 
therefore not the right quantity to measure the size of a single nucleon. Unfortunately for this simple 
trial action we have no alternative. For our other trial actions the y-particle describes purely self 
energies and therefore the distance of a particle from its y-particle is an unambiguous contribution 
to the cloud radius. We will discuss this point further when we present our results. As in the one 
particle case the cloud radius for this simple trial action is always larger than the rms (core) radius 
approaching each other in the strong coupling limit. For N ^ 1 the cloud radius approaches the limit 

^cioud -^'^^P^ (3.38) 

which is equal to the interparticle distance. This is to be expected since the definition given in (3.37) 
includes the exchanges which probe the interparticle distance. We remark that what we call here 
cloud radius is in agreement for N = 1 with the polaron radius as defined in [3]. 

All the corresponding formulae for the above quantities in the case of the 3- and 4-parameter trial 
actions are given in Appendix B. 
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4 Results and Discussion 



As it was mentioned in the introduction we are interested in the qualitative description of the nuclear 
many body system mainly due to the lack of antisymmetry. We are thus addressing questions like 
what is the relative importance of self energies and vertex corrections in a saturating bound system? 
How do these change as we go to a strongly bound system? What is the mean number of mesons and 
the amount of dressing in a saturating system and how are they modified as the binding is increased? 
The purpose is not to reproduce the binding energy of real nuclei nor to describe accurately the 
Deuteron where one would need a much more sophisticated calculation. Therefore our approach in 
this work is not to search the whole parameter space in order to fix the scalar and vector couplings 
and form factors to reproduce the binding energies of real nuclei. Instead we consider four sets of 
these parameters that we regard representative with respect to the posed questions. The values of the 
couplings gs and gv entering the Walecka model, and of the cutoffs As and Ay are listed in table 1 
and we denote these sets by Gi,G2,G3 and G4. Gq in table 1 is a set of parameters taken from 
ref. [15] and resulted from a fit to nuclear properties in a mean field approximation. Go has no form 
factors and leads to divergent loop integrals. We therefore introduce a scalar and vector form factor 
each with a cut off. Following ref. [3] we fixed the bare mass of the nucleon M by requiring that the 
mass of a single "nucleon" be equal to the physical nucleon mass i.e. we have set the effective mass 
Mefr(l) = 938.9 MeV. As was pointed in [3] increasing the value of the form factors decreases the 
bare mass to an extent that our nonrelativistic approximation can no longer be justified. Therefore 
the values of the cut offs are restricted for a given set of couplings. For large couplings the cut offs 
become close to the meson masses. Decreasing the couplings allows an increase in the values of the 
cut offs. With these considerations we compromised with the sets listed in table 1. There we also 
give the values for the bare mass and other single particle observables. For N = 1 trial action 1, 2 
and 4 coincide and give the results of table 1 whereas trial action 3 produces very similar values ^. 
The two first sets with the stronger couplings reduce the bare mass by more than 600 MeV and are 
at the limit of applicability of the non-relativistic nucleon dynamics. G3 and G4 on the other hand 
produce reasonable bare masses. 

The form of the static two-body potential corresponding to these sets of parameters is presented 
in fig.l. Including the form factors for the couplings of Go softens the core and reduces the depth of 
the potential. As we can see in fig. 1 the maximum attraction is the same for the two sets Gi and G2. 
The repulsive core progressively decreases as we go from Gi to G4. The ratio of maximum attraction 
to maximum repulsion increases from 0.03 for Gi to 0.15 for G3 to the large value of 1.8 for G4. 

Having fixed the bare mass we determine the variational parameters by minimizing the energy 
E]\f of the N-particle state ^. The dependence on N of the variational parameters obtained in the 
minimization of the N-particle energies is displayed in fig. 3 for all trial actions and for the parameter 

^ Using trial action 3 decreases the values of the bare mass by 12 MeV in the case of parameter sets Gi and G2 and 
by 0.4 MeV for G3 and G4. 

^The minimization was performed with the CERN subroutine MINUIT while the expressions for the energy and the 
effective mass were evaluated numerically with typically 64 Gauss points. 
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sets G2 and G'4 using lines for the former and symbols for the latter ^ . The parameters W£, and V£, of 
the purely diagonal terms in Sp are much larger than v and w and therefore W£, and V£, are scaled 
by 0.1 in order to display them on the same figure. For action 2 we have set w = W£, and since v > w 
all parameters are scaled by 0.1. For set G2 the dependence of all the variational parameters on N 
is very weak and although the strength C of the exchange term is small, including this term in the 
action lowers the energy. For the set G4 the N dependence of w and v is stronger and the strength 
C becomes larger, in agreement with the fact that the exchange terms dominate the behaviour of 
a strongly bound system. Having determined the variational parameters we proceed to discuss the 
various observables. 

In table 2 we list the energies obtained from the minimization using the four trial actions for the 
parameter set G'2. The biggest improvement comes from including a coupling of each nucleon to a 
fictitious particle modeling the self energy which trial action 1 lacks The variational bound is lower 
by a significant amount by using trial action 2 indicating the importance of having different strengths 
C and C'd for the direct and exchanged terms even though the retardation is the same. Allowing in 
addition a different retardation for the exchange terms as is done for trial action 3 and 4 lowers the 
energy by a smaller amount. The fact that trial action 3 produces the best bound, also compared to 
trial action 4, indicates the importance of the extra diagonal term of strength the same as that of the 
exchange terms a feature shared by the effective action, ^eff. 

Since the one-particle energies are the same for trial actions 1,2 and 4 the variational principle 
applies on the level of their binding energies per particle defined as 



The results for the binding energy per particle as a function of N are summarized in fig. 4 for the 
four different sets of parameters and trial actions. For the set Gi the system exhibits an almost 
saturating behaviour but, for these values of the couplings and with no antisymmetry, with an order 
of magnitude larger BE]\f as compared to the nuclear system. For the rest of the parameter sets the 
slope becomes increasing larger with G4 as an extreme case where BE]\f grows rapidly. The main 
observation here is that whereas trial action 1 gives no binding for the sets with the larger repulsion 
(Gri, G'2) it yields progressively better results as we go to G3 and G4. Similarly results from trial 
actions 3 and 4 become progressively more similar as we go to G3 and G4, whereas results using trial 
action 2 always coincide with those from trial action 4. Therefore we conclude that for a saturating 
system it is very important to construct a good trial action whereas for a strongly bound system 
this becomes less crucial. The fact that trial action 2 yields results the same as 4 indicates that it 
is more important to allow a different strength for the diagonal and non-diagonal terms rather than 
a different retardation. Having compared the various trial actions we proceed to discuss the results 
obtained using trial action 3 which is superior to the rest. 

It is interesting to look separately at the diagonal and non-diagonal contributions to the total 

■^The behaviour of the variational parameters for the sets Gi and G3 is in overall similar to that of G2. 
* Allowing for two different exponentials (4-parameters) in trial action f , which corresponds to coupling all nucleons 
to two fictitious Z— particles, leads to a minor improvement of the variational bound for N > 1. 



BEn 



En - NEi 
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binding energy. These are shown in fig. 5 for all four sets of parameters using trial action 3. For Gi 
and G2 we observe a large cancelation between the diagonal and exchanged contributions. As the 
repulsive core is decreased the slope of the exchanged contribution increases changing from negative to 
positive at a decreasing value of N . The diagonal contribution remains constant and therefore the N 
dependence of BE]\f is determined by the behaviour of the exchanged part which becomes dominant 
for the set G4. 

We can further break up the total diagonal and non-diagonal contributions to those coming from 
A^* = i < ^eff — Sf > and those coming from the energy, Ep^ corresponding to the trial action. 
These separate contributions to the energy and binding energy per nucleon are plotted in fig. 6 for 
the set Gi and G4 using again trial action 3. We notice that in both cases the diagonal part of Ep is 
small and constant and therefore the contribution of this term to BEn is just its single particle value. 
In the case of Gi the rest of the terms grow linearly with N and therefore give a constant contribution 
of the same order of magnitude to BEn ^ . Going from the set Gi to G2 we begin to see a non-linear 
dependence of the non-diagonal part of A^* which leads to the slope of BE^. This observation is 
independent of the trial action and is due to the increase of the binding in our model. When we go 
to G'4 the non-diagonal A^* piece has a significant slope and becomes the dominant contribution. For 
this last set all trial actions produce results in much closer resemblance to those displayed in fig. 6. 

The effective mass defined by (3.22) is a measure of the kinetic mass of the iV-particle state which 
in this non-relativistic limit has no direct connection to the energy. In analogy to the binding energy 
we may define a mass defect as 

If the mesons were real particles in the N-particle state then we could define a mass defect in terms 
of their masses 

Afeff(iV) - {NM+ < Ns> ms+ < Ny > my) 



AMn = ^ ' ^ ■ f ^ (4.41) 



which, compared to AMjv, is shifted by AMi. These are given in table 1 for our four parameter sets. 
Shifts of the order of the nucleon mass point again to the breakdown of the non-relativistic treatment 
of nucleons. One can also define an excess meson mass by 

msANs + myANv 

AMs,v = — ^^^^ (4.42) 

where 

AN, =< N,(N) > -N < N,(l) > i = S,V (4.43) 

are the meson excess in the N-particle system. Whereas the effective mass was calculated in first order 
the meson numbers are calculated to zeroth order. The results are plotted in fig. 7 for all parameter 
sets. The mass defect increases more rapidly with N than the excess meson mass. The vector excess 
per particle, AiVy/iV, is almost zero for the sets Gi and G2 rising rapidly for G4. As expected there 
are, on the average, much more scalar mesons present. 



^The results obtained with trial action 2 resemble very closely those of trial action 3 whereas for trial action 4 the 
values of AS are much smaller leaving as major cancelling contributions the diagonal and non-diagonal parts of Ep- 
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Finally we estimate various sizes of the N-particle state, all of which are calculated to zeroth order 
using the intrinsic wave function derived from the underlying model Hamiltonian. In fig. 8 we show 
the interparticle distance which settles to about 2 fm for Gi and, as expected, decreases to about 0.5 
fm for G4. We generalize the definition of rem given for trial action 1 in section 3 to apply to the 
other actions with the extra fictitious particle attached to each nucleon. We take 



where the nucleon feels the presence of all other nucleons through the Z— particle. For trial action 3 
this follows the interparticle distance expect for small N where the system is unstable with respect 
to two-particle break-up. 

How one measures the extension of the meson cloud around one bare nucleon is not very well 
defined since the Z— particle is shared among all the nucleons. A lower bound is given by considering 
just the y— particle i.e. the expectation value of (xj- — yi)^. For trial action 4 this represents the 
only self energy contribution and it was the main motivation for introducing this trial action. As it 
can be seen this cloud radius is constant and small. An upper bound is given by summing both the 
contribution from the y— and the Z— particles i.e < (Z — Xj)^ -|- (xj- — yi)^ >. The latter follows 
closely the interparticle distance. 

5 Summary and Conclusions 

We have presented a variational calculation of the N-particle bound state where dynamical meson 
exchange is summed up to all orders. We have for the first time included all vertex corrections and 
self energies in an N particle system albeit in an approximate way. While the scalar and vector 
mesons are treated fully relativistically the N constituents are treated non-relativistically. The main 
advantage of this method is that we are able to obtain analytic expressions for the observables of the 
N-body bound state and therefore study the behaviour as a function of N . 

The crucial step in this work as in all variational approaches is the construction of a good trial 
action. For this we used the original idea by Feynman of replacing mesonic degrees of freedom 
by fictitious particles. This study revealed that it is of great importance to separate direct and 
exchange contributions giving them different strengths. This turns out to be much more important 
than distinguishing vector and scalar contributions. Studying trial action 4 we learn that although 
the strength of the purely exchange term in this action is small compared to that of the direct it 
significantly lowers the energy bound obtained with trial action 1 which makes no distinction between 
the self-energies and exchanges. 

Using the Walecka model as our underlying theory requires the introduction of form factors. 
For the typical couplings of this model one has to tune the form factors to avoid large self energy 
contributions. We found that for a particular set of parameters we approximated a saturating system 
where the exchanged terms contribute a constant to the binding energy per particle cancelling to a 
large extent with the diagonal part. For these large cancelations occurring in a nearly saturating 




N 



(4.44) 
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system it is important to construct a good trial action. This becomes less of an issue for a strongly 
bound system where the exchange contributions dominate. Another important feature that emerges 
from this calculation is that the number of excess mesons per constituent in a saturating system is 
rather small being practically zero as far as vector mesons are concerned and reaching about two for 
scalar mesons when N ~ 100. The radius of the dressed constituents defined solely in terms of the self 
energy comes out rather small and independent of the amount of binding. In contrast the interparticle 
distance is more than five times larger for the saturating system becoming almost equal to the cloud 
radius in the large binding situation. 

Although the above features can not be taken as a quantitative description of the real nuclear 
case due to the lack of antisymmetrisation, nevertheless they are useful in that they reveal important 
difference between a saturating and a strongly bound system. In order to approach the nuclear system, 
in addition to the issue of antisymmetry, one has to treat the it — N interaction which is strongly spin 
and isospin dependent. Treating the more realistic Chew-Low model [16] using coherent spin states 
[17] gives rise to a complex effective action [18] requiring a better understanding before variational 
results of similar accuracy as those obtained in the polaron problem can be produced. There has been 
some work [19] recently on the treatment of spin in the particle representation but it remains an open 
question whether we can apply those ideas within the present framework. 

Acknowledgement: We thank R. Rosenf elder for many illuminating discussions and for a careful 
reading of the manuscript. 
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A Determination of the static effective potential 



We derive here the static two-body potential from the effective action ^eff. Setting Fi(k) = 1 in eq. 
(2.8), we can perform the momentum integration in terms of the modified Bessel functions [20] leading 
to the following expression for the effective action 



« j'=i 

N 



1=1 

fi N 



ms 



— / dt' \^ rriy / da Kq I mv\/a'^\^i(t) — ^j(t' 



ay 



+ \t-t 



l\2 



N 

ay \ - 
2^ ^ 



1 



(A.l) 



where constant terms coming from tadpole diagrams are omitted since they are irrelevant for the two- 
body potential. Due to the exponential suppression 'm_T = t — t' the main contribution in (A.l) comes 
from trajectories with Xj(i) ~ Xj(i') ~ Xj(r) where T = In this approximation and neglecting 

the term proportional to the transverse velocities we can perform the r integration [20]. We get the 
following result for the non-diagonal part of the action 
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Using the relation K^if2i^) 
well known result 



e ^, the static potential between a pair of nucleons reduces to the 



V{r) = ay- 



-rmy 



-TVfls 



as- 



T — X -j' 



(A.3) 



r r 

a superposition of a Yukawa attraction due to scalar meson exchange and a Yukawa repulsion due to 
vector meson exchange. 
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B Calculation of the relevant observables for a general quadratic 
action 



We give here the details for the calculation of the expectation values needed to extract physical 
observables. They are applicable to any action quadratic in the particle trajectories. 



B.l Energy and effective mass 

The main quantity which has to be calculated is the mean value < e*-(^'(^)-^^(<^)) >Sp with respect 
to the trial action Sp. For this purpose we apply the stationary phase approximation to the action 

rf3 N 



S' = S- [ dt J2 W-Mi 

•^0 a = l 



where the source fa(^) is given by 

Ut) = ik(6,Mt - r) - 6,,6(t - a)) . 
Using the trajectories ii(t) which make S' stationary one easily obtains 

'l ff^ 



< e 



jk-(X,(T)-Xj((7)) 



>Sp = exp 



1 ^ 

^ a = l 



:b.i] 



:b.3) 



In this way the problem reduces to finding Xj-. The required mean value follows by inserting the 
resulting expression in the above formula. To simplify the calculations we work in Fourier space using 
the following decomposition for the trajectories of the particles 



Mt) = ^^a? + E cos(27rn-) + ^ sin(27rn-; 

^ n=l P n=l P 



:b.4) 



where, when two indices are used, subscripts denote particle labels and superscripts the Fourier 
component. In this representation we get, for a general quadratic action, the following set of equations 
which the Fourier coefficients a°,a" and must satisfy in order to make the action S' stationary 
(we suppress the tildes in the following) 
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Due to translation invariance we are free to choose Aq = 0. The coefficients A„,A„,Ao,Ao are in 
general rational functions of n and their exact form depends on the form of the trial action we use. 
For trial action 3 and 4 we obtain respectively 
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The expressions for trial actions 1 and 2 can be obtained from those of trial action 3 by setting cd = 
and w = W£) respectively. 

The system of linear equations given in eq. (B.5) has the solution 



Xn 
m 

vr>, _j_ 

W° 



A X" 



A 



iVA. 



Ao 



fB.lO) 



Substituting in (B.4) these stationary values we obtain from (B.3) 
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where 
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valid in the limit /3 ^ 00. We have define the sums Si and S2 by 
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In order to perform these sums we have to first find the roots in the numerator of these expressions 
and decompose the corresponding rational functions of n in partial fractions. From this point on 
we will consider quadratic actions with any sum of exponential retardation factors. In this case the 
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numerator A„ or A„ — iVA„ can always be written as a real polynomial P in ra^ of degree p and the 
sums for the corresponding Sj in the above expressions are given by 
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where —p'ji are complex numbers denoting the roots of P and Cj^i are complex coefficients resulting 
from the partial fraction decomposition. To determine the contribution to the sum Sj coming from a 
conjugate pair of roots we have to be able to calculate the following types of sums 
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Since we are interested in the limit /3 ^ oo of the above expressions we can replace the sums by 
integrals and use the integration formulas of [20] to obtain Ti and T2. Having determined the 
sums and consequently the mean value < e*-(^'(^)-^^(<^)) 

>Sp we can repeat the above 

procedure using a suitable source term (see [3] for details) to evaluate the mean value < x_L^j(r) • 



jk-(X,(T)-Xj((7)) 



>Sp- The rest of the expectation values needed in (3.21) can be obtained 
from (B.3) by suitable differentiation with respect to k. It is then straightforward to obtain also the 
effective mass using the definition of (3.22). The results for the energy and the effective mass are 
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where we used the notation 

if i = J 
if i 7^ J 

The functions G{t) and H{t) are given by the second derivative of Si(r) and S2(t). The constant 
Co is given by eq.(3.28) and is independent of the variational parameters. The expressions for Eq and 
< Sp >Sp in (B.15) depend on the exact form of the trial action and they are given explicitly below 
for trial action 3. The quantities Ri(k,T) are defined in eq.(2.9). 
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In general the calculation of ^2i't) can only be accomplished numerically. For trial action 

3 however the polynomial P reduces to a biquadratic form and we can perform the calculations 
analytically to obtain for the energy of the iV-particle system 

1 pf3 poo r / 2 2 
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and the mass 
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The quantities l-i^{T) and I^DiT) are given by 
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and the functions G(t) and -ff (t) by 

G(T) 

a ■ ' c- 

The quantities i7, 77, are determined in terms of the variational parameters through 
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The ground state energy Ep for trial action 3 reads 



rB.22) 
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which vanishes in the free case. 

We can again recover the results for trial action 1 by taking cp, Q and for trial action 2 by 
taking c ^ in the above formulae. For trial action 4 these expressions for the energy and effective 
are evaluated numerically. 



B.2 Radii 

The mean inter-particle distance fjv is given in general by 

1 ^ 3 

< rl >= — ^ E < - >'= ^ MO) (B.24) 

The calculation of the cloud and rms radius is in general more complicated and the resulting expres- 
sions depend on the form of the trial action. In general one has to resort to a numerical evaluation. 
We describe here the calculation of these quantities in the case of trial action 3. To define the rms or 
center of mass radius, rem, we consider the subsystem of the ith nucleon, the Z— and the yi particles. 
We take the distance of the nucleon from the center of mass of this subsystem to be the rms radius 
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The cloud radius should give a measure of the extension of the meson cloud around the bare nucleon 
and therefore one should sum up the all the self energy contributions. For the single particle system 
this is given by < (Z — x)^ + (y — x)^ >i/2 when a y-particle is present. For more than one particle 
the Z-particle does no longer correspond to a purely self energy contribution but models also the 
exchange of mesons between different nucleons. Since the y-particles describe purely the self energies 
we may take 
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as a lower bound to the cloud radius. The quantity 
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greatly overestimates the dressing since it probes via meson exchange also the interparticle separation. 
In trial action 4 the contribution of the Z— particle to the self energy is turned off by hand by simply 
not including the corresponding diagonal term. Although for this action the cloud radius is given 
entirely by < >i/2 we do not have the wave function since this action is not described by a model 
Lagrangian. Using the wave function of trial action 3 with the variational parameters derived for 4 is 
the approximation which we adopt. 

To calculate the expectation values in the above equations we find the wave function of the system 
described through the Langrangian (3.17). The easiest way to do this is to go over to normal mode 
coordinates. For this we use the scaled variables x' = \/Mxi, y' = ^Jrn^yi and Z' = yjMzTj and we 
express the potential energy part of (3.15) in terms of these variables. The corresponding potential 
energy matrix for trial action 3 is given by 
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We then determine numerically the matrix A which brings V into a diagonal form. In terms of 
this matrix the coordinates x',y', Z' are given by 
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where q^, i = 1, ..,2N + 1 are the normal mode coordinates of the system. In these coordinates the 
wave function of the system is given by 
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where are the eigenvalues of the matrix V . Due to translation invariance one of these eigenvalues 
corresponding to the center of mass vanishes. Without loss of generality we take the vanishing 
eigenvalue to be ei. Then the expectation values in eq. (B.25-B.27) with respect to the wave function 
of (B.30) can be calculated without difficulty leading to the following expressions 
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In the case of trial action 1 the above expressions can be evaluated analytically leading to the formulae 
given in section 3. 
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Figure captions 

Figure 1 : The static potential of eq.(2.14) for five different values of the scalar and vector couplings and 
cutoffs. The dotted line corresponds to Go (infinite cut offs and as = 6.52 and ay = fO.84), the 
long dash-dotted to d (as = 6.52, A5 = 0.54 GeV and ay = f 0.84, Ay = 0.92 GeV), the long 
dashed to G2 (as = 3.2, A5 = f.O GeV and ay = 5.f2,Ay = f.4 GeV), the short dash-dotted 
to G3 {as = 2.4, A5 = f.O GeV and ay = 4.0, Ay = f.O GeV), and the short dashed to G4 
(as = 3.4, A5 = f.O GeV and ay = 3.57 Ay = f.2 GeV). 

Figure 2 : The system of harmonic oscillators which describes the trial action of eq.(3.f8). 

Figure 3 : The variational parameters of the trial actions f to 4 are shown in (i) to (iv) respectively. We 
use lines to show the parameters for the set G2 and symbols for the set G4. The purely diagonal 
variational parameters are scaled in both cases by a factor of O.f. They are shown by the 
dotted line ( retardation Wd/10) and by the short dashed line (vd/10), whereas the long dashed 
(retardation w) and dash-dotted (v) lines show variational parameters for the exchange terms. 
The corresponding diagonal parameters are shown by the triangles and crosses and the exchange 
by the squares and x's. In (iv) we do not plot the variational parameters for the set G4 since 
their behaviour is the same as for the set G4 in (iii). 

Figure 4 : The binding energy per particle in GeV versus the number of particles. For the parameter set 
Gi trial action f produces no binding, trial action 3 and 4 are shown with the dash-dotted and 
solid lines respectively. For set G2 trial action f produces no binding, trial action 3 and 4 are 
shown with the dashed and dotted lines respectively. The crosses and squares correspond to 
trial action f and 3 (same as trial action 4) for the set G3 and the stars and triangles to trial 
action f and 3 (same as trial action 4) for the set G4. The results for trial action 2 coincide in 
all cases with those from action 4. 

Figure 5 : The binding energy per nucleon as a function of N for trial action 3 is shown by the dashed 
line. The squares and crosses show the diagonal and non-diagonal contributions, 
(a) Parameter set Gi, (b) Parameter set G2, (c) Parameter set G3, (d) Parameter set G4. 

Figure 6 : (a) The energy in GeV is shown versus the number of particles by the dashed line for the parame- 
ter set Gi using trial action 3. The squares and triangles denote the diagonal contributions from 
A^* and the trial energy respectively and the x's and crosses the corresponding non-diagonal 
parts, (b) The same as (a) but for the parameter set G4. 

(c) The same as (a) but for the binding energy per particle. 

(d) The same as (c) but for the binding energy per particle. 

Figure 7 : The effective mass defect (triangles) and the excess meson mass (crosses) per particle in GeV 
as well as the scalar (dash-dotted line) and vector (dotted line) meson excess per particle are 
shown using trial action 3 for the different parameter sets. The (a)-(d) correspond to the same 
set of values as those denoted in figure 5. 
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Figure 8 : The interparticle distance is shown with the dotted line and the center of mass radius with 
the squares both using trial action 3. The cloud radius is shown with the dashed line and it 
is obtained using the parameters of trial action 4. The notation (a)-(d) is the same as that in 
figure 5. 
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Table 1: We give the values of the scalar and vector couplings and cut ofFs, gs, gv, and Ay 
respectively and denote a given set of these values by the symbol given in the first column. The rest 
of the columns give single particle properties which coincide when using trial actions 1,2 and 4. The 
last six columns give the values of the bare mass and the mass defect in GeV, the average number of 
scalar and vector mesons and the rms and cloud radius in fm. 





as 


av 


A5 


Ay 


Mo 


AMi 


<Ns> 


< Nv > 




^cloud 








(GeV) 


(GeV) 


(GeV) 


(GeV) 






(fm) 


(fm) 


Go 


6.52 


10.84 


CXD 


CXD 
















6.52 


10.84 


0.54 


0.921 


0.181 


-0.888 


0.21 


1.96 


0.24 


0.28 


G2 


3.2 


5.12 


1.0 


1.4 


0.272 


-1.193 


0.21 


2.23 


0.18 


0.24 


Gs 


2.4 


4.0 


1.0 


1.0 


0.578 


-0.267 


0.16 


0.69 


0.12 


0.36 


G4 


3.4 


3.57 


1.0 


1.2 


0.517 


-0.469 


0.22 


0.99 


0.13 


0.30 



Table 2: The values of the total energy Ej^' and the binding energy BEj^' for the parameter set G2 
using the trial action j,(j = 1, ..,4) are given for various particle numbers N . 



N 


^(1) 


^(2) 


^(3) 


^(4) 












(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


1 


-1.790 


-1.790 


-1.838 


-1.790 


0.0 


0.0 


0.0 


0.0 


2 


-2.482 


-3.656 


-3.741 


-3.677 


-0.549 


0.038 


0.032 


0.049 


4 


-4.423 


-7.311 


-7.482 


-7.325 


-0.684 


0.038 


0.032 


0.042 


6 


-6.431 


-10.975 


-11.222 


-10.987 


-0.718 


0.040 


0.032 


0.042 


8 


-8.485 


-14.649 


-14.963 


-14.646 


-0.730 


0.041 


0.032 


0.041 


10 


-10.549 


-18.338 


-18.727 


-18.335 


-0.735 


0.044 


0.035 


0.044 


20 


-21.273 


-37.088 


-37.872 


-37.130 


-0.726 


0.065 


0.055 


0.067 


30 


-32.545 


-56.359 


-57.541 


-56.402 


-0.705 


0.089 


0.080 


0.091 


40 


-44.352 


-76.157 


-77.739 


-76.192 


-0.681 


0.114 


0.105 


0.115 


50 


-56.691 


-96.483 


-98.467 


-96.535 


-0.656 


0.140 


0.131 


0.141 


100 


-126.333 


-206.048 


-210.045 


-206.188 


-0.526 


0.271 


0.262 


0.271 
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Figure 2 
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